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Abstract 

We present a measurement of the composition of B meson inclusive semileptonic decays using 
9.4 fb _1 of e + e~ data taken with the CLEO detector at the T(45) resonance. In addition to 
measuring the charged lepton kinematics, the neutrino four- vector is inferred using the hermiticity 
of the detector. We perform a maximum likelihood fit over the full three-dimensional differential 
decay distribution for the fractional contributions from the B — > X c iv processes with X c = D, D*, 
D**, and nonresonant X c , and the process B — > X u tv. From the fit results we extract the first 
and second moments of the M\ and q 2 distributions with minimum lepton-energy requirements 
of 1.0 GeV and 1.5 GeV. We find (M\ - M 2 D ) = (0.456 ± 0.014 ± 0.045 ± 0.109) GeV 2 /c 4 with a 
minimum lepton energy of 1.0 GeV and (M\ - M* D ) = (0.293 ± 0.012 ± 0.033 ± 0.048) GeV 2 /c 4 
with minimum lepton energy of 1.5 GeV. The uncertainties are from statistics, detector systematic 
effects, and model dependence, respectively. As a test of the HQET and OPE calculations, the 
results for the M\ moment as a function of the minimum lepton energy requirement are compared 
to the predictions. 

PACS numbers: 13.20.He, 12.15.Hh, 12.39.Hg 
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I. INTRODUCTION 



Inclusive semileptonic B meson decay can be used to measure the CKM parameters 
involved in the decay and to measure nonperturbative hadronic properties of the B meson. 
Heavy quark effective theory (HQET) combined with the operator product expansion (OPE) 
provides a framework in which many inclusive B decay properties can be calculated In 
particular, moments of the differential decay rates of a variety of processes are related 
to nonperturbative parameters that also appear in the calculation of the total decay rates. 
Measurements of these moments can therefore be used to refine calculations of the B — ► X u lv 
and B — > X c £u decay rates and the extraction of the CKM parameters \V u b\ and \V c b\ from 
measurements of the respective branching fractions. In this paper, we present measurements 
of the first and second moments of the M\ and q 2 kinematic variables in the decay B — ► X c lv^ 
where M\ and q 2 are the squares of the invariant masses of the hadronic and leptonic parts 
of the final state, respectively. 

In the HQET and OPE framework, the inclusive B decay matrix elements are expanded 
in powers of Aqcd/^b- For each order in the expansion new nonperturbative parameters 
arise: at order Aqcd/Mb, there is A; at order Aq CD /M^, there are Ai and A2; and at 
order Aq CD /M^, there are pi, p 2 , T 2 , T 3 , and T±. The nonperturbative parameter A 
relates the b quark mass to the B meson mass in the limit of infinite b quark mass. The 
Ai and A2 parameters are related to the kinetic energy of the b quark inside the B meson 
and the chromomagnetic moment of the b quark inside the B meson. The parameter A2 
is directly related to the mass splitting between the B* and B mesons. These parameters 
describe properties of the B meson and are not specific to the decay mode being studied. 
For example, the same parameters appear in expansions of the moments of the lepton energy 
and X c mass distributions in B — > X c iv decays and in expansions of the moments of the 
photon-energy spectrum in B — ► X s ^ decays. 

These calculations do not predict the long-distance effects that govern the formation of 
hadrons. They are therefore only applicable when a sufficiently large region of phase space is 
included in an observable that the hadronization effects are negligible. The differential decay 
rates themselves do not satisfy this condition. Instead, moments of the differential decay 
rates are measured and compared to the HQET-OPE predictions. The assumption made in 
applying quark-level calculations to hadron-level processes is known as quark-hadron duality. 
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We report an analysis of B — > X£v decays in which both the charged lepton and the 
neutrino are reconstructed. We use a maximum likelihood fit to the full three-dimensional 
kinematic distribution to extract the exclusive branching fractions for B — > D£u, B — > D*£u, 
B — > D**£u, nonresonant B — > X c £u, and i? — > X u £v. The descriptions of the kinematic 
distributions of these components used in the fit are derived from theoretical calculations or 
models. The extracted branching fractions are very sensitive to these models. However, the 
description of the inclusive differential decay rate constructed from these models and the 
extracted branching fractions is less model dependent, because the branching fractions have 
been adjusted by the fit to replicate the measured inclusive distribution. The (M| — M D ), 
((M| — (M 2 .)) 2 ), (g 2 ), and ((q 2 — (q 2 )) 2 } moments of the B — > X c £v differential decay rate 
where Mr> is spin averaged D meson mass, Md — (Mp + 3 x Md*)/4, are calculated from 
this description and are similarly less model dependent than the branching fractions. 

The moments measured in this paper can be used to determine the HQET-OPE pa- 
rameters. At present, the knowledge of the A and Ai parameters is sufficient so that the 
uncertainty in the extraction of \V c b\ from the semileptonic B meson decay rate is due to the 
uncertainty in the contributions from the third-order terms and second-order a s corrections. 
The main goal of further moments measurements is to over-constrain the determination of A 
and Ai, and thereby test the quark-hadron duality assumption. These moments could also 
be used to constrain the third-order terms to further improve the precision of the extrac- 
tion of \V c b\. The dependence of the moments on the minimum lepton energy used in the 
measurement provides an additional test of the theoretical assumptions and consistency. 

II. DETECTOR AND DATA SET 

The data used in this analysis were taken with two configurations of the CLEO detector, 
CLEO II and CLEO II. V. An integrated luminosity of 9.4 fb _1 was accumulated on the T(AS) 
resonance (E cm ps 10.58 GeV), and an additional 4.5 fb _1 was accumulated at 60 MeV below 
the T(AS) resonance, where there is no BB production. Both detector configurations covered 
95% of the An solid angle with drift chambers and a cesium iodide calorimeter. Particle 
identification was provided by muon chambers with measurements made at material depths 
of 3, 5, and 7 hadronic interaction lengths, a time-of-flight system and specific ionization 
(dE/dx) measured in the drift chamber. In the CLEO II configuration, there were three 
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concentric drift chambers filled with a mixture of argon and ethane. In the CLEO II. V 
detector, the innermost tracking chamber was replaced with a three-layer silicon detector 
and the main drift chamber gas was changed to a mixture of helium and propane. The 
CLEO II and II. V detectors are described in more detail in references y and y. 



III. EVENT SELECTION AND RECONSTRUCTION 

Events are selected to have an identified electron or muon with momentum greater than 
1 GeV/c and a well- reconstructed neutrino. Additional criteria are used to suppress back- 
ground events from the e + e~ — > qq continuum under the T(4S') resonance, where q is u, d, 
s, or c. 

The identified leptons are required to fall within the barrel region of the detector 
( | cos 6*| < 0.71, where 9 is the angle between the lepton momentum and the beam axis). 
Electrons are identified with a likelihood-based discriminator which combines dE/dx, time- 
of-flight, and the ratio of the energy deposited in the calorimeter to the momentum of the 
associated charged track (E/p). Muons are identified by their penetration into the muon 
chambers. For momenta between 1.0 and 1.5 GeV/c, muon candidates are required to pene- 
trate at least 3 interaction lengths and above 1.5 GeV/c, candidates are required to penetrate 
at least 5 interaction lengths. The lepton-identification efficiencies are calculated by embed- 
ding raw data from reconstructed leptons in radiative QED events into hadronic events. 
The rate at which pions and kaons fake leptons is measured by reconstructing Kg — > 7r + 7r~, 
D° — > K~7i + , and D° — > K + n~ using only kinematics and then checking the daughter 
particle lepton-identification information. 

Neutrinos are reconstructed by subtracting the sum of the four-momenta of all observed 
tracks and showers not associated with tracks, £>ot> served , fr° m the four- momentum of the 
e + e _ initial state, p^+ e ~, which is nearly at rest in the laboratory: 

Vv Pe + e~ ^observed - 

The errors made in this assumption are due to particles lost through inefficiency or limited 
acceptance, fake tracks and showers, and other undetected particles such as K Q L mesons, 
neutrons, or additional neutrinos. Several requirements are made to select events in which 
these effects are reduced and the neutrino four-momentum resolution is correspondingly 
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improved. 

Because extra neutrinos are correlated with extra charged leptons, events with an identi- 
fied lepton in addition to the signal lepton are rejected. The primary source of fake tracks is 
charged particles that do not have sufficient transverse momentum to reach the calorimeter 
and therefore curl in the tracking chambers, returning to the beam axis. The drift chamber 
hits produced by such a particle after its initial outbound trajectory may be reconstructed 
as additional tracks. Criteria have been developed to identify such errors and make a best 
estimate of the actual charged particles in the event. Events for which the net charge of 
the all tracks selected by these criteria is not zero are removed, reducing the effect of lost 
or fake tracks. Showers in the calorimeter that are matched to tracks in the drift chamber 
are not used, to avoid double-counting their energy. A neural network algorithm provides 
further rejection of secondary hadronic showers associated with showers that are matched 
with tracks. 

A final neutrino reconstruction quality requirement is that the mass of the reconstructed 
neutrino must be small. The ratio of the reconstructed neutrino invariant mass squared 
to twice the reconstructed neutrino energy is required to satisfy \M*/2E V \ < 0.35 GeV/c 4 . 
This quantity is approximately proportional to the energy of a lost or fake particle. After 
this cut, the reconstructed neutrino's energy is assigned to be the magnitude of the missing 
momentum, because the momentum is not dependent on the particle identification of the 
tracks and so has a better resolution than the direct energy measurement. The resulting 
neutrino energy resolution has a narrow core with a full width at half maximum height of 
approximately 120 MeV and a broad tail of over estimation of the neutrino energy which 
extends up to 1.5 GeV. 

Continuum events are suppressed by a combination of event-shape and -orientation cri- 
teria. These exploit the tendency of continuum events to be jet-like and aligned with the 
beam axis, in contrast with BB events which are more spherical and randomly oriented in 
the detector. The ratio of second to the zeroth Fox-Wolfram moment jj], i?2 = H2/H0, of 
the energy flow in the event is required to be less than 0.4. In addition, a neural network 
is used to combine R2, the angle between the lepton and the event thrust axis, the angle 
between the lepton momentum and the beam axis, and the fraction of the total energy lying 
in nine separate cones around the lepton direction, which cover the full Atc solid angle. The 
i?2 cut is more than 99% and 95% efficient for B — > X c iv and B — > X u £u, respectively, while 
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removing 60% of the continuum events. The neural-net cut removes an additional 73% of 
the continuum background, while keeping 92% and 94% of B — > X c iv and B — > X u £u, 
respectively. 

After all cuts we observe 121851 events from the data sample collected on the T(4S) 
resonance. The overall efficiency varies from 1.5% for B — ► X c lv nonresonant to 4.2% for 
B -> XJv. 

IV. KINEMATIC VARIABLES 

The differential decay rate of inclusive semileptonic B meson decays can be described in 
terms of three independent kinematic variables, which can be chosen to be the squares of 
the masses of the hadronic and leptonic parts of the final state [M\ and q 2 ) and the cosine 
of the helicity angle of the virtual W (cos Owe)- The helicity angle is defined as the angle 
between the lepton momentum in the virtual-W^ frame and the virtual- W momentum in the 
B meson frame. 

Because we do not reconstruct the hadronic part of the final state, M\ must be inferred 
through kinematics: 

M 2 X = M 2 B + q 2 - 2E hcam {E £ + E v ) + 2\p B \ \p t + p v \ cos 6 B . £u , 

where q = pe+p u is the momentum of the lepton system, p B is the momentum of the B, and 
Qb iv is the angle between them. Since the B mesons are the daughters of an T(45 l ) produced 
at rest, the magnitude of the B momentum is known and small (\pb\ ~ SOOMeV), however 
its direction is unmeasured. The last term the M\ formula depends on the B momentum 
direction, and is small, unmeasured, and neglected in this analysis. Because of the neglected 
term, the M\ resolution depends on In addition, because of the unknown B momentum 
direction, the B meson frame in the definition of cos 6wt is replaced by the lab frame in the 
definition of the reconstructed quantity. 

V. COMPOSITION EXTRACTION 

The full three-dimensional differential decay rate distribution as a function of the recon- 
structed quantities q 2 , M\, and cos#vw is fitted for the contributions from semileptonic B 
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decay and backgrounds. The q 2 variable is replaced by q 2 /{En + E v ) 2 for fitting purposes. 
This has the effect of varying the q 2 bin size as a function of Eg, + E v . The B decay modes 
are B -> Dip, B -> £>*.&/, B -> D**£v, B -> X c £i/ nonresonant, and B -> XJz/. The back- 
grounds are classified as secondary leptons, continuum leptons, or fake leptons. A secondary 
lepton is a real lepton in a BB event whose parent is not a B meson. A continuum lepton 
is a real lepton in a continuum e + e~ — > gg event. A fake lepton is a non-leptonic track from 
either a BB or a continuum event which is identified as a lepton. 

We perform a binned maximum-likelihood fit where component histograms are con- 
structed from weighted Monte Carlo or data events. The fit uses electrons and muons 
simultaneously, with a separate set of histograms for each. The likelihood is implemented 
to take into account the histogram statistics using the method described in reference jsj. 

The B —>■ Xiv modes, secondary leptons and real leptons from the continuum are modeled 
with events from a GEANT 71 simulation of the CLEO detector that are reconstructed in 



the same manner as data events. The B — > Dtv and B — > D*lv modes are simulated with an 
HQET-based model using the PDG |3| averages of measurements of the form factors rescaled 
to have the curvature term set to 50% of its theoretically predicted value (see Section IVlDJ) . 
The B — > D**£i> and B — ► X u iv modes are simulated using form factors from the ISGW2 
model P|. The X c nonresonant modes are simulated with the Goity and Roberts model 10 1 
in which the D and D* contributions are excluded. 

The fake leptons are modeled with data events where a track is selected to be treated as 
a lepton. The events are unfolded bin by bin to extract the ir and K contributions, which 
are then multiplied by the measured fake rates. This models fake leptons from both BB 
and continuum. This method also provides an absolute normalization for the fake-lepton 
contribution to the data sample. The real leptons from the continuum are modeled with a 
Monte Carlo simulation which has been tuned to replicate the appropriate charm spectra; 
charm is the source of most leptons from continuum. The models of both continuum and 
fake leptons have been validated and constrained by comparisons with the 4.5 fb" 1 of off- 
resonance data. The secondary leptons are modeled with CLEO's generic BB Monte Carlo 
which has also been tuned to replicate measured charm spectra and semileptonic charm- 
decay measurements. The measured branching fractions for B meson decays to charmed final 
states do not sum to the theoretical prediction for the inclusive rate [llj . This discrepancy 
referred to as the charm counting problem is most likely due to missing or mismeasured 
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modes in the sum. The branching fractions in the BB Monte Carlo simulation are therefore 
tuned to saturate the theoretically predicted level of charm production. 

Final-state radiation (FSR) can play an important role in semileptonic decays. This 
is particularly important for events with an electron in the final state, because the small 
electron mass enhances the effect. For the B — ► X r iv modes, the Monte Carlo events 
sifted with GEANT indude the effects of radiation ^ the PHOTOS pa ckag e Q to 
generate radiated photons. This simulated both the physics of radiation and the detector 
response to the photons. For the B — > X u iv modes and the e + e~ — > qq continuum Monte 
Carlo events, the events were simulated without FSR and an algorithm to apply the FSR 
after the detector simulation is used. This algorithm generates photons and calculates the 
effect on the lepton four-vectors in the same way as the PHOTOS package. The change 
in the lepton momenta is then applied to the reconstructed leptons. The effect of losing a 
photon is simulated using a random number to apply the photon efficiency extracted from 
the GEANT simulation. If the photon is rejected, its four-vector is added to the neutrino 
four-vector simulating the effect of the additional lost particle on the reconstructed neutrino 
kinematics. 

The normalization of the continuum lepton component is determined from the data taken 
below BB threshold. The normalization of the fake leptons is determined from the measured 
fake rates and the measured track spectra. The contributions of these two backgrounds are 
therefore not allowed to vary in the fit, while those of the secondary leptons and all of the 
B — > Xiv modes are. A summary of the processes contributing to the selected sample, the 
fraction of the sample each contributes, and the models used to describe them in the fit is 
shown in Table HJ The contribution fractions are either determined by the fit or externally 
constrained as described above. 

Projections of the Monte Carlo simulations of reconstructed quantities q 2 , and 
cos Owe for the various B — > Xiv modes are shown in Figure H These one-dimensional pro- 
jections illustrate some of the discriminating power available to the full three-dimensional 
fit. Projections of the data and fit result, shown in Figure 121 are compatible with each other 
within the estimated size of the systematic uncertainties. 

From the fit results the branching fractions are determined using the efficiency as pre- 
dicted by the Monte Carlo simulation, the data yield, and the number of BB pairs in the 
data sample. The resulting branching fractions are shown in Table ITT1 
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TABLE I: Composition of the data sample and summary of the models used in the fit. 



Mode 



Fraction of Data Sample Model 



B -» Dlv 
B -» D*£u 
B -» D**lv 

Nonresonant B — > X c ^i/ 
S -» 

Secondary Leptons 

Fake Leptons 

Electron 



Ef, < 1.5 GeV 
£ M > 1.5 GeV 



Continuum Leptons 



0.118 
0.476 
0.084 
0.033 
0.016 
0.050 

0.132 
0.002 
0.100 
0.030 

0.089 



HQET [5J 
HQET [5] 
ISGW2 [9] 

Goity and Roberts [10] 
ISGW2 [9] 

CLEO B decay model & measurements 

of semileptonic charm hadron decay 
Data & measured lepton fake rates 



JETSET [12J 



VI. SYSTEMATIC UNCERTAINTIES 

The method of neutrino reconstruction adds a large amount of kinematic information to 
each event. However, it also adds significant potential for systematic errors. The Monte 
Carlo simulation of resolution on the neutrino kinematics is affected by the modeling of 
the signal, the other B in the event, and the detector response. The GEANT Monte Carlo 
simulation does not perfectly reproduce the track and shower efficiencies and fake rates, 
nor are B decays well enough understood that the inclusive particle distributions are well 
known. For this analysis we employ a reweighting method in order to quantify the effects of 
these uncertainties on our results. For example, to study the effect of the tracking efficiency 
uncertainty, the Monte Carlo events in which tracks are lost are given a higher or lower 
weight in constructing the component histograms used in the fit. 

The shifts of the nominal result due to variations of the detector performance, the model- 
ing of inclusive B and D meson decay, and radiative corrections are summed in quadrature 
to get the total detector systematic uncertainty. The shifts for each variation of the signal 
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FIG. 1: Monte Carlo simulation of the distributions of the reconstructed quantities q 2 , M^, and 
cos#vK£ f° r the various B — ► Xtv modes. The curves are normalized to have unit area to facilitate 
comparison of the shapes. Note that due to imperfect resolution M\ can be less than zero. 

B model are similarly summed in quadrature to get the total model-dependence systematic. 
The larger contributions to the systematic uncertainties are summarized in Table IIHI 

A. The Modeling of Inclusive B and D Decay 

The probability that the full set of final-state particles in an event will be found depends 
on the number of particles, their momenta and their type. The model of the B and D 
decay physics used in the Monte Carlo is tuned to reproduce a wide variety of inclusive and 
exclusive measurements. We have identified a number of inclusive properties to which the 
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*3100304-012 




Nonresonant 
B— X c £v 

□ Backgrounds 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
Neutrino Energy (GeV) 



FIG. 2: Projections of the fit results in the reconstructed quantities. On the first row are q 2 
and M\] on the second are cos#iy^ and Eg, and on the third row is E v . The histograms are the 
Monte Carlo simulation and the points with error bars are the data. The step in the lepton-energy 
distribution is due to the looser muon identification used below 1.5 GeV/c for which there is a 
higher fake rate and higher efficiency. 
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TABLE II: Branching fractions results. The errors on the entries in the table are the statistical, 
detector systematic, and model dependence uncertainties, respectively. 



Mode Branching Fraction (x 10 2 ) 

B — > Dtu 1.92 ± 0.08 ± 0.19 ± 0.74 

B -» D*iv 6.37 ± 0.06 ± 0.65 ± 0.86 

B -» D**tv 1.51 ± 0.07 ± 0.30 ± 0.52 
Nonresonant B -> Jf c lz/ 0.70 ± 0.07 ± 0.25 ± 0.62 

B -» 0.12 ± 0.01 ± 0.03 ± 0.23 

Sum 10.61 ± 0.29 ± 1.09 ± 0.75 



neutrino resolution and the efficiency of the event selection are particularly sensitive. These 
are the number of K\ mesons, baryons, and extra neutrinos in the final state, and the total 
charged particle and photon multiplicities. 

The extra neutrinos in the events come predominantly from semileptonic decays of the 
other B meson in the event, and from secondary charm-decays, B — > c — > x£v. These 
are both suppressed by the second-lepton veto, but because of the energy threshold for the 
lepton identification, there is a significant unvetoed contribution. The branching fraction 
B(B — ► c — > xiv) calculated from the measured charm production and charm semileptonic 
branching fractions is 9.6 ±0.9% jfij. However, the branching fractions for charm production 
in B meson decay are not consistent with the theoretically predicted rate The inclusive 
B decay model used in the simulation, which is consistent with the theoretical prediction, 
gives a higher rate of 10.7%. Based on these numbers the systematic uncertainty assigned 
to this rate is ±1.1%. 

The number of baryons in B meson decay is determined from measurements of the branch- 
ing fraction B(B — ► (p or p)X) = 8.0 ±0.4% Because of the uncertainty in the exclusive 
composition of the process, the number of B — ■> baryons events in the simulation is varied 
by ±20%. 

The number of K\ mesons is inferred from a measurement of the inclusive number of Kg 
mesons in B decay. The discrepancy between this measurement and the value used in the 
simulation is 7%. The number of K\ mesons is therefore conservatively varied by ±10%. 
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TABLE III: Contributions to the systematic uncertainties of the branching fraction measurements. 
The entries are given as a percentage of the central value. The entries separated by a slash indicate 
the effect of raising and lowering the varied quantity, respectively. The uncertainties quoted for 
the B — > X c tv nonresonant mass distribution and B —> X u iv models represent the range covered 
by the set of models studied. 











Nonresonant 




Inclusive 


Variation 


B — ► DiV 


B — > D*iv 


B — ► D**£u 


B — ► X c £v 


B — » X u £u 


B — ► Xlv 


Statistical 


3.9 


0.9 


4.9 


10.1 


7.0 


2.7 


Detector 


10.1 


10.2 


20.0 


35.2 


25.5 


10.3 


Model Dep. 


38.7 


13.5 


34.2 


88.7 


201.5 


7.0 


Lcpton Fake Rate 


0.8 / -0.6 


0.4 / -0.3 


2.7 / -1.8 


6.9 / -2.8 


12.6 / -6 .3 


0.8 / -0.2 


Continuum ±10% 


-3.6 / 2.6 


0.4 / -0.5 


-2.9 / 6.7 


1.7 / -6.3 


-4.1 / 4.6 


-0.8 / 0.7 


B(b — > c -» t) ±10% 


1.9 / -2.0 


2.6 / -2.5 


6.6 / -6.3 


-0.8 / 0.7 


-6.1 / 6.4 


2.7 / -2.6 


B(b -* baryons) ±20% 


1.7 / -0.9 


3.9 / -3.0 


4.9 / -3.9 


1.7 / -1.0 


-0.1 / 0.9 


3.4 / -2.6 


# K l 


2.4 / -2.4 


2.6 / -2.5 


7.0 / -6.7 


-1.3 / 1.2 


-7.2 / 7.7 


2.8 / -2.7 


Track Efficiency 


-4.8 / 4.9 


-5.8 / 6.0 


-9.6 / 10.8 


0.9 / -1.7 


-1.7 / 1.7 


-5.6 / 5.9 


# Fake Tracks 


1.9 / -1.7 


2.3 / -2.0 


2.6 / -2.3 


-3.9 / 3.5 


0.3 / -0.2 


1.9 / -1.6 


Shower Efficiency 


-2.2 / 2.5 


-1.4 / 1.7 


-3.9 / 4.8 


1.8 / -1.4 


-3.8 / 3.7 


-1.7 / 2.1 


# Fake Showers 


-3.1 / 2.7 


-0.5 / 0.2 


4.5 / -3.9 


-28.2 / 27.5 


-14.7 / 15.3 


-2.2 / 2.0 


Track Multiplicity 


2.5 


2.0 


3.1 


2.2 


-3.0 


2.2 


n r\~\\T£>v l\/T 1 1 1 "t" i t~\ 1 1 r*i ^~~\r 
OIl*J Wcl 1VI U.1L1 jJULl Ly 


2.2 


3.4 


3.8 


8.9 


-3.4 


3.5 


Final-State Radiation 


-0.1 


-3.3 


-0.2 


9.0 


-7.2 


-1.5 


Lepton Efficiency 


4.3 


2.0 


4.1 


5.1 


-3.7 


2.8 


B — > Dlv p Param. 


6.4 / -5.2 


-2.6 / 1.8 


4.6 / -2.0 


-5.4 / 2.3 


0.1 / 0.2 


-0.1 / 0.0 


B — > ZWf c Param. 


1.4 / -1.1 


-0.5 / 0.3 


0.8 / -0.4 


-0.5 / 0.2 


0.1 / -0.0 


0.1 / -0.1 


B -> D*&/ p Param. 


32.4 / -32.1 


-11.4 / 11.0 


17.9 / -15.2 


-18.0 / 17.0 


6.3 / -4.8 


0.4 / -0.3 


_B — > D*&/ Param. 


-1.1 / 0.8 


0.8 / -0.7 


-1.0 / 0.6 


2.9 / -2.0 


1.1 / -1.0 


0.3 / -0.3 


B -> D*iv R1/R2 1st Eig. Vec. 


-13.3 / 14.9 


3.5 / -3.7 


-10.0 / 10.1 


16.7 / -15.3 


-4.5 / 5.7 


-0.7 / 1.0 


B -► D*&< R1/R2 2nd Eig. Vec. 


9.9 / -10.4 


-4.0 / 4.2 


3.2 / -3.4 


0.9 / -0.7 


-3.3 / 3.5 


-0.1 / 0.2 


_B -> HQET Model 


-4.0 


0.5 


2.3 


0.1 


0.2 


-0.1 


B -> D**£u w Slope 


-2.7 / 3.1 


0.9 / -0.8 


-2.7 / 0.6 


13.6 / -9.8 


-0.5 / 0.2 


0.5 / -0.5 


B — » X c <?f Nonres ui Slope 


-0.0 / -0.1 


0.1 / -0.0 


-6.3 / 3.6 


17.4 / -9.9 


-0.2 / 0.1 


0.3 / -0.2 


.B — » X c fo Nonres Mass 


1.0/-6.2 


1.0/-2.9 


25.6/-4.9 


81.8/-13.5 


2.5/0.2 


6.9/-2.7 


B -> X u lu Model 


3.0/-3.0 


0.2/-1.5 


0.8/-2.3 


0.9/-6.7 


201.2/-19.4 


0.5/-0.2 



To assess the uncertainty due to the charged particle and photon multiplicities of in- 
clusive B decay, simulated events are reweighted to correctly reproduce track and shower 
multiplicities observed in the data sample. The full shift is used as an estimate of the 
uncertainty. 
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B. The Detector Response 



The effect of the detector response enters through the lepton-identification efficiency 
and fake rate, and the neutrino energy resolution and efficiency The lepton-identification 
efficiency is varied between the efficiencies determined by the embedding procedure described 
above, and the Monte Carlo simulation. The shift is larger than the statistical errors on the 
embedding measurement. This conservatively estimates the effect of a systematic shift, and 
shows lepton efficiency to be a very small contribution to the uncertainty. 

The uncertainty due to lepton fake rates is determined by using Gaussian distributed 
random numbers to vary the fake rates within their experimental errors. These variations 
are then propagated through the fake lepton model for a maximum of ten trials. The shift 
of the results due to an increase or decrease of the total rate is also included. 

The neutrino-momentum resolution and event-selection efficiency are strongly affected 
by how well the tracks and showers observed correspond to the actual number of charged 
particles and photons produced in the event. This correspondence is affected by the track 
and shower efficiency and fake rate. 

The uncertainty on the charged particle efficiency is determined using the constraints 
of charge conservation in 1-prong versus 3-prong r pair events to infer when a track has 
been lost. Monte Carlo simulated tracks embedded in hadronic data events are used to 
study the effect of the event environment on the tracking efficiency. The uncertainty on 
the photon efficiency is determined by varying the inputs into the Monte Carlo simulation 
and by comparing the shower shape distributions observed in data with the Monte Carlo 
simulation. Both of these Monte Carlo studies are described in detail in reference 

Tracks identified by the reconstruction software that do not correspond to actual charged 
particles produced in the collisions are referred to as fake tracks. The uncertainty of the 
Monte Carlo prediction for the number of fake tracks is estimated using 1-prong versus 
3-prong r pair events with an extra track and the total charge distribution of the event 
sample. 

The largest detector-related uncertainty is the number of reconstructed showers that are 
not due to photons. The main cause of these showers is secondary hadronic interactions of 
the particles produced in the primary showers of charged hadrons. The reliability of the 
Monte Carlo simulation of these showers has been studied in 77 — > K^Kg events and in 



16 



t — » 7T 7rV T events. In both cases, the expected number of photons in the detector is 
well defined for the specified mode. The number of observed showers is then compared to 
the number of showers predicted by a Monte Carlo simulation. These measurements are 
imprecise because of the uncertainty in the contributions of other modes to the selected 
sample. Based on these comparisons, a variation of the number of fake photons by ±10% is 
used to assign an uncertainty. 



C. Radiative Corrections 



PHOTOS implements an algorithm based on a splitting function that applies the same 
physics at 0(a) as the prescription of Atwood and Marciano Q|. PHOTOS also modifies 
the kinematics of the decay to force the conservation of momentum in addition to energy. 
The PHOTOS algorithm implements an approximation that ignores the internal structure of 
the hadronic system. Richter-Was j.16] has made a comparison of the PHOTOS and Atwood- 
Marciano prescriptions with an exact order-a calculation of the radiative corrections to the 
B~ — > D°£~V differential decay rate and has found agreement at the 20% and 30% level, 
respectively. Because we must extrapolate to the other B — > X c iv modes, we make a 
conservative estimate that the PHOTOS calculation can be trusted to ±50%. 



D. Signal Mode Model Dependence 

The models for all the B — > Xtv components were varied to assess the model dependence 
of the measured branching fractions. The B — > Dlv and B — > D*iv modes were varied 
within the range of the errors on the form-factor measurements. The curvature of the form 
facto, has not been and is nsna„ y confined b y theoretic*, predictions Q. 

Because the data have an excess above the model in the q 2 region between 5.0 and 8.0 GeV 2 , 
the curvature was set to 50% of its predicted value and varied by ±50% of the prediction. 
The ISGW2 model for the B — > D**tv form factors [^| was replaced by a model inspired 
by HQET calculations 18]. The slope of the B — > D**£u and nonresonant B — > X c tv form 
factors in the q 2 dimension was also varied. 

The dominant model dependence is from the hadronic mass distribution of the B X c iv 
nonresonant mode. This is conservatively reweighted with a series of Gaussians restricted 
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to the kinematically allowed region. The means of the Gaussians are allowed to range from 
+ M n to 3.5 GeV/c 2 , with variances ranging from 0.25 GeV 2 /c 4 to 1.25 GeV 2 /c 4 . 

The B — > X u lv simulation is varied from an all-nonresonant model to the nominal ISGW2 
model Q], with a hybrid of the two in between. The all-nonresonant model differential decay 
rate corresponds to the pre 
B — > Xg'y spectral function 

The maximum deviation of the nonresonant B — > X c iv mass Gaussian variations is 
added in quadrature with the deviation of the other model variations to get the total model 
dependence. 



iction of HQET combined with CLEO s measurement of the 
19, Q- 



VII. CALCULATION OF MOMENTS 



The branching fraction of the individual components, B m , combined with the physics 
models used in the fit form a description of the inclusive differential decay rate. It is from 
this description that the moments results are calculated. A moment (M) with a cut C is 
expressed as 

{M) c = y c B ' ( ) 

where c m is the fraction of the decay rate for mode m in the region defined by the cut, 

f ^C(x)d x 
J ^jf- dx 

and m m is the moment of the mode m in that region 



JM(x)^C(x)dx 
mm - J^C(x)dx ■ [6) 

The quantities m m and c m depend only on the model. The measured branching fractions, B m , 
depend on the model, the detector simulation, and the data. Since the branching fractions 
are measured using the inclusive differential decay rate, when combined with the models 
used in the fit they give a good description of the true differential decay rate. Mismodeling 
of a contribution may cause the branching fraction to be mismeasured, but the shape will 
still be well described. For instance, the main separation of the D and D* modes is due to 
the q 2 distribution. If the q 2 slope of either of these modes is mismodeled, the relative rates 
of the modes will be affected, but the model of the q 2 distribution and its moments will be 
only weakly affected. 
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A. Radiative Corrections 



Radiative corrections play an important role in the measurement of the M\ distribution. 
The reconstructed M\ is defined to be the system recoiling against the charged lepton and 
the neutrino. If a photon is radiated by the lepton in the event, it will be included in this 
definition of the recoil system, 

Mi = (^-^-^) 2 = (K + ^) 2 - (4) 

The goal of this analysis is to measure the mass-squared moment of the recoiling hadronic 
system p\, not its combination with the radiated photon (ph+P-y) 2 - To correct for this effect, 
the data are fit using fully simulated GEANT Monte Carlo events in which the PHOTOS 
package [13J has been used to generate radiated photons. The moments are calculated from 
the fit results and the models without radiative corrections of the hadronic mass distributions 
for each mode and thus do not have a shift due to the radiative corrections. 

The application of the radiative corrections is further complicated by the fact that the 
generated photons are low energy and often lost. When the photon is lost it can cause the 
event to fail the missing- mass cut. If the event does pass the missing-mass cut, the recon- 
structed neutrino will be biased toward high energy, pushing the reconstructed M\ toward 
the true hadronic mass without the photon. If neglected, this would increase the measured 
(M x — M D ) moment with a 1.0 GeV lepton-energy cut by 0.082 GeV 2 /c 4 , before detector 
effects are included. This is reduced to 0.037 GeV 2 /c 4 after detector effects. The variation 
used to assign the systematic uncertainty due to the radiative corrections is discussed in 
Section ESI 

B. Results 

The resulting (M x - M 2 D ), ((M x - (M x )) 2 ), (q 2 ), and ((q 2 - (q 2 )) 2 ) moments of the 
B — > X c iv differential decay rate with their uncertainties are presented in Table IIV1 The 
corresponding correlation coefficients are shown in Table |3 These correlations include 
statistical, systematic, and model-dependence uncertainties. Because the moments measured 
with a 1.0 GeV and a 1.5 GeV lepton-energy cut are highly correlated, it is also useful to 
consider the difference between the (M x — M D ) moment with the two different cuts. We 
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find 



(M 2 x ) Ee>1 . 0GeV - (Mx) E t >i.5GeV = (0.163 ± 0.014 ± 0.036 ± 0.064) GeV 2 /c 4 

and a correlation coefficient of this value with the (M x — M D )Eg>i.5GeV moment of 0.486. 
The (M x - M D ) moments SIS db function of the lepton energy are shown in Table IVII The 
contributions of the individual systematic uncertainties for the moments results with a 1.0 
GeV lepton-energy cut are shown in Table IVTT1 

TABLE IV: Moments results with E e > 1.0 GeV and E e > 1.5 GeV lepton-energy cuts. The 
errors on the entries in the table are the statistical, detector systematic, and model dependence 
uncertainties, respectively. 



Moment 
(M X -M 2 D ) (GeV 2 /c 4 ) 
{(Ml - (Mi)) 2 ) (GeV 4 /c 8 ) 
(q 2 ) (GeV 2 ) 
((g 2 - (g 2 )) 2 ) (GeV 4 ) 



E t > 1.0 GeV 
0.456 ± 0.014 ± 0.045 ± 0.109 
1.266 ± 0.065 ± 0.222 ± 0.631 
4.892 ± 0.015 ± 0.094 ± 0.100 
2.852 ± 0.002 ± 0.003 ± 0.047 



E e > 1.5 GeV 
0.293 ± 0.012 ± 0.033 ± 0.048 
0.629 ± 0.031 ± 0.088 ± 0.113 
5.287 ± 0.020 ± 0.073 ± 0.095 
2.879 ± 0.006 ± 0.007 ± 0.049 



TABLE V: Correlation coefficients of the moments measurements presented in Table llVl 



Moment 


Cut (GeV) 


Correlation Coefficients 


(M|- 


-M 2 D ) 


E t > 1.0 GeV 


1.000 


0.910 


0.970 


0.881 


-0.795 -0.651 


-0.034 


-0.122 


(M x - 


-M 2 D ) 


E e > 1.5 GeV 




1.000 


0.824 


0.856 


-0.814 -0.784 


-0.103 


-0.179 


((M 2 X 


~ (M x )?) 


E t > 1.0 GeV 






1.000 


0.884 


-0.683 -0.523 


0.054 


-0.036 


((M 2 X 


~ (Ml)) 2 ) 


E t > 1.5 GeV 








1.000 


-0.606 -0.531 


0.116 


0.052 


(q 2 ) 




Ei > 1.0 GeV 










1.000 0.925 


0.301 


0.352 






E t > 1.5 GeV 










1.000 


0.406 


0.442 




(q 2 )) 2 ) 


Et > 1.0 GeV 












1.000 


0.979 


((q 2 - 


(q 2 )) 2 ) 


Et > 1.5 GeV 














1.000 
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TABLE VI: {M%- — M D ) versus the lepton-energy cut. The errors on the entries in the table are 
the statistical, detector systematic, and model-dependence uncertainties, respectively. 

Cut (GeV) (M\-M 2 D ) (GeV 2 /c 4 ) 

E t > 1.0 0.456 ± 0.014 ± 0.045 ± 0.109 

E e > 1.1 0.422 ± 0.014 ± 0.031 ± 0.084 

Et > 1.2 0.393 ± 0.013 ± 0.027 ± 0.069 

Et > 1.3 0.364 ± 0.013 ± 0.030 ± 0.054 

E t > 1.4 0.332 ± 0.012 ± 0.027 ± 0.055 

E t > 1.5 0.293 ± 0.012 ± 0.033 ± 0.048 



VIII. CONCLUSION 

We have measured the first and second moments of the M| and q 2 distributions in the in- 
clusive process B — > X c £u. Results are presented with minimum lepton-energy requirements 



of 1.0 and 1.5 GeV. We also present the (M\ — M D ) moment as a function of the lepton- 
energy cut. The M\ moments with the 1.5 GeV lepton-energy requirement are in good 
agreement with CLEO's previously reported M\ moments j^jj, obtained with the CLEO-II 
subset of the data used here. The results reported here, in addition to using more data, 
handle final-state radiation more carefully. These results supersede the previously published 
results. 

The M\ and q 2 moments reported here, moments of the lepton-energy spectrum previ- 
ously reported 22], moments of the lepton energy spectrum over a broader energy range 
23] , and moments of the photon-energy spectrum in the radiative penguin decay B — > X s ~f 
20j can all be interpreted in the context of the HQET-OPE framework. The measurements 
should collectively provide a good determination of the HQET-OPE nonperturbative pa- 
rameters A and Ai, provide constraints on the C(Aq CD /M^) nonperturbative parameters 
Pi, P2-, %.i %, and T4, and provide a test of the quark-hadron duality assumption. An 



interpretation of this body o 
the errors, is in preparation 
presented in this paper. 



CLEO data, taking proper account of the correlations among 
Here we give two examples of interpretation of the results 
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TABLE VII: Contributions to systematic uncertainties of the moments measurements with a 1.0 
GeV lepton-energy cut. The uncertainties due to the B — » X c lv nonresonant mass distribution 
and B — ► X u lv models represent the range covered the set of models studied. 





(M x - M D ) 


(i M x ~ \ M x)> ) 




// 2 / 2\\2\ 

(w -(<?})> 


Variation 


(GeV /C J 


(GeV /e°J 


1 r> T (2 \ 

(GeV ) 


(in — 3 /~i "\7"4\ 

(10 ° X GeV ) 


Lepton Fake Rate 


0.014 / -0.006 


0.047 / -0.019 


-0.011 / 0.006 


-0.840 / 0.362 


DELCO b — > c — ► I shape 


-0.008 


-0.017 


0.010 


0.421 


DELCO +1ct b — > c — » £ shape 


-0.011 


-0.021 


0.015 


0.574 


DELCO -Ict — > c — » £ shape 


-0.001 


-0.005 


-0.001 


0.036 


Continuum Norm ±10% 


0.001 / 0.001 


0.006 / -0.029 


0.015 / -0.020 


0.149 / -0.198 


p, Fakes, Ei < 1.5 GeV, ±10% 


-0.011 / 0.011 


-0.068 / 0.070 


0.005 / -0.003 


0.803 / -0.793 


p Fakes, Et > 1.5 GeV, ±10% 


-0.000 / 0.000 


0.009 / -0.009 


0.001 / -0.001 


-0.044 / 0.049 


B(b — * c — > ±10% 


0.004 / -0.003 


-0.015 / 0.015 


-0.006 / 0.006 


-0.068 / 0.048 


S(6 — > baryons) ±20% 


0.002 / -0.002 


-0.011 / 0.011 


0.003 / -0.003 


0.096 / -0.104 


# ^ 


n nno / n nn o 

0.003 / -0.003 


n n 1 o / n m n 

-0.018 / 0.019 


n nn o 1 n nn o 

-0.008 / 0.008 


n c\c 1 / nnoo 

-0.061 / 0.033 


Irack h/tnciency 


n nno / n nm 

0.002 / -0.001 


n n A n / c\ r\ a ~\ 

0.040 / -0.041 


n nno / n nn a 

0.003 / -0.004 


n noo / n nrr 1 ? 

-0.288 / 0.257 


# rake 1 racks 


n /~inn / n nn o 

-0.009 / 0.008 


n n A i / n noo 

-0.041 / 0.038 


n nric / n nn a 

0.005 / -0.004 


r\ a o / n rr nc 

0.548 / -0.506 


Shower Efficiency 


n nn 1 / n nnn 

0.001 / 0.000 


nmo / n n i c 

0.018 / -0.016 


n n n c I n nnT 

0.006 / -0.007 


n n / n m o 

-0.069 / -0.013 


# Fake Showers 


n nnn / n no-i 

-0.030 / 0.031 


-0.181 / 0.172 


n n-i o / n n 1 >i 

0.013 / -0.014 


n "i nn / n i >i 

2.129 / -2.164 


lorce Irk Multiplicity 


0.002 


0.003 


-0.004 


-0.120 


lorce bnwr Multiplicity 


0.011 


0.037 


-0.003 


-0.534 


Final-State Radiation 


0.021 


0.088 


-0.024 


-1.567 


Lepton Efficiency 


0.006 


0.022 


-0.012 


-0.499 


Total Detector 


0.045 


0.222 


0.094 


3.130 


B — > Diu p Param. 


-0.006 / 0.005 


-0.016 / 0.004 


0.022 / -0.010 


-7.112 / 5.057 


B — > D£u C£) Param. 


-0.001 / 0.001 


-0.000 / -0.000 


0.000 / 0.001 


5.850 / -4.784 


B — > D*iv p Param. 


-0.020 / 0.025 


-0.045 / 0.046 


0.041 / -0.059 


26.712 / -27.140 


B — > D*iv ca x Param. 


0.003 / -0.002 


0.015 / -0.010 


-0.010 / 0.006 


26.477 / -26.716 


B -> D*lv R1/R2 1st Eig. Vec. 


0.019 / -0.018 


0.081 / -0.070 


-0.047 / 0.044 


-19.883 / 19.427 


B -> R1/R2 2nd Eig. Vec. 


-0.000 / 0.000 


0.028 / -0.028 


-0.006 / 0.005 


-1.072 / 0.859 


B — * D**lu HQET Model 


0.008 


0.007 


-0.014 


10.384 


B -> «; Slope 


0.014 / -0.014 


0.077 / -0.058 


-0.013 / 0.015 


12.122 / -11.871 


B — > X c fo( NonRes Slope 


0.008 / -0.005 


0.050 / -0.031 


-0.006 / 0.003 


2.837 / -1.456 


B — > X c fe( Nonres Mass 


0.102/-0.078 


0.615/-0.543 


-0.058/0.029 


5.928/-4.284 


B -> Model 


0.009/-0.012 


0.018/-0.039 


-0.004/0.005 


-0.339/0.580 


Total Model 


0.109 


.631 


0.100 


0.047 



In Figure El we present a comparison of our results for the (My — M D ) moment as a 
function of the lepton-energy cut with HQET-OPE predictions [25]. The A and Ai param- 
eters are constrained by the first photon-energy moment of the B — > X s ^ process and 
the (JW| — M D ) moment with a 1.5 GeV lepton-energy cut from this analysis. The theory 
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FIG. 3: (M| - M D ) Versus Lepton-Energy Cut. The "CLEO 2003" data points are from the work 
presented here and the "CLEO 2001" data point is from Ref. j21| . 

bands shown in the figure reflect the experimental uncertainties on the two constraints, the 
variation of the third-order HQET parameters by the scale (0.5 GeV) 3 , and variation of the 



size of the higher order QCD corrections 



25| . The theoretical calculation and the results of 
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this measurement appear to agree. The more detailed analysis, in preparation 
the correlations of the measurements will provide a more stringent test of the prediction 
of the HQET-OPE theory. A violation of quark-hadron duality could manifest itself as a 
discrepancy in this calculation. 

In Figure 01 we show the bands in A-Ai space defined by the (Mf - M D ), the {{M 2 X - 
(Mjr)) 2 ), and the (Mx)E t >i.oGeV ~ (^x)^>i.5GeV moment measurements. The widths of 
the bands reflect the experimental uncertainties on the measured quantity, the variation of 
the third-order HQET parameters by the scale (0.5 GeV) 3 , and variation of the size of the 
higher order QCD correct ions [3]. As shown Table W\ there are strong correlations among 
the errors, and hence a simple band plot only gives a qualitative indication of the values of 
A and Aj. A precise determination, with errors, awaits the full analysis in preparation [24] . 
The extracted values of these parameters, when combined with precision measurements of 
the B XA branch haetio, Q m d the B mes on MeW Q, pe rml t the extract™ 
of \Vcb\ with reduced theoretical uncertainties. 

During the final preparation of this paper, we learned of a preprint from the BaBar 
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FIG. 4: Constraints on the nonperturbative parameters A and Ai due to the {M\ — M D ) and 
((Mj^ — (M^-)) 2 ) moment measurements with minimum lepton-energy requirement of 1.0 GeV and 
the {Mx)E e >i.oGeV — (Mx)E e >i.5GeV moment measurement (left). Constraints on the nonpertur- 
bative parameters A and Ai due to the (M\ — M 2 D ) and ((Mj^ — (M^)) 2 ) moment measurements 
with minimum lepton-energy requirement of 1.5 GeV (right). 

collaboration reporting new meaurements of the first four hadronic mass moments . The 
second and fourth hadronic mass moments reported by BaBar correspond to the first and 
second hadronic mass-squared moments reported here and are consistent within the quoted 
uncertainties. The BaBar measurements are based on an approximately nine times larger 
T(AS) sample and use a significantly different analysis technique. Their technique results in a 
larger statistical uncertainty and a smaller systematic uncertainty. The combined statistical 
and systematic uncertainties of the BaBar measurements range from half the uncertainty 
quoted here to similar uncertainty for the various comparable measurements. 
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